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Abstract. Ultracompact stellar models with a two-zone uniform density equation of state 
are considered. They are shown to provide neat examples of optical geometries exhibiting 
double necks, implying that the gravitational wave potential has a double well structure. 
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1 Introduction 

The static spherically symmetric (SSS) relativistic stellar models have been the subject 
of extensive studies in the past (see for example for a review of exact perfect fluid 
SSS solutions of the Einstein equations up to 1994). De Felice ||] first discussed the 
possibility of trapping of massless particles in the interior of compact stellar objects. 
Abramowicz and coworkers || considered the optical geometry of uniform density 
models. They found that the optical geometry (for geodesies in the equatorial plane) 
of the stellar interior was identical to a standard 2-sphere. They also showed that the 
optical geometry of the uniform density models with R < 3M had a neck (at R = 3M) 
and a corresponding bulge in the interior. The bulge is associated with stable spatially 
closed orbits of massless (or ultrarelativistic) particles and waves propagating at the 
speed of light (in the geometric optics approximation). This provided a nice visual 
illustration of the phenomenon of trapping of waves and matter moving at relativistic 
velocities inside compact objects. Moreover, it allowed a simple derivation of the 
eigenfrequencies of the w-modes of trapped gravitational waves. 

Despite the relative simplicity of the field equations for static spherically symmet- 
ric perfect fluid models, the full picture of the structure of the gravitational field of 
such objects is only now beginning to emerge. Following the powerful Nilsson-Uggla 
formulation of the SSS field equations a novel feature which was recently uncov- 
ered is the existence of optical geometries with multiple necks || . One might naively 
expect that such behavior would only occur for rather strange or unphysical equations 
of state. However, as shown in [^J, multiple necks do in fact appear for the innocent 
looking Zel'dovich stiff matter equation of state with an added constant p = p — p s . 
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We are using a subscript "s" to denote evaluation at zero pressure (i.e. at the stellar 
surface). In fact if n is the number of the necks, then n — > oo for stellar models with 
that equation of state in the limit of infinite central density, p c — > oo, where the sub- 
script "c" is used for evaluation at the center of the stellar object. Actually, most of 
these models correspond to unstable equilibria and so would not be realizable as stable 
stellar objects. However, unstable equlibria can play a role as intermediate states in 
gravitational collapse situations Q. 

The Zel'dovich equation of state has the advantage of being causal in the sense 
that the adiabatic speed of sound equals the speed of light, u SOU nd : = y dp/ dp = 1. 
However, from the point of view of visualizing the multiple neck optical geometry, it is 
not so well suited because of the rather small amplitude of the necks for that equation 
of state. Therefore, to obtain better visual displays we have used double layer uniform 
density models to exhibit double neck optical geometries. Such models were used by 
Lindblom as prototypes to discuss phase transitions in relativistic stellar objects. 



2 Some multiple neck optical geometries 

The metric for a SSS spacetime can be written in the Schwarzschild radial gauge as 

ds 2 = -e 2v dt 2 + e 2X dr 2 + r 2 (d9 2 + sin 2 9d(/) 2 ), (1) 

where v and A are functions of r. When studying null geodesies, it is often convenient 
to introduce the optical geometry of the spacetime, defined by 

ds 2 = e- 2v ds 2 = -dt 2 + e 2(x -^dr 2 + e - 2 V# 2 . (2) 

Since this metric is conformally related to the physical metric, the spacetime and its 
optical geometry have the same null geodesies. Due to the symmetries of the spacetime 
we may, without loss of generality, consider only the spacelike slice given by t — and 
8 = 7r/2. Introducing the tortoise radial variable defined by dr* = e x ~ v dr : the metric 
on this slice may be written 

dl 2 = drl+r 2 d<p 2 , (3) 

where f — e~ u r. In order to gain intuitive understanding of this geometry it is useful 
to embed it in 3-dimcnsional Euclidean space. Standard embedding techniques (c.f. 
e.g. ||) yield the differential equation for the embedded surface 



r = e-"r, ti = e~> 'yj e 2A - (1 - rv') 2 , (4) 

where f and h are the radius and height coordinates in a standard cylindrical coordi- 
nate system and a prime denotes differentiation with respect to r. 

As shown by Chandrasekhar and Ferrari the axial modes of nonradial metric 
perturbations of SSS perfect fluids are (to first order) not coupled to a perturbation 
of the matter. Such metric perturbations may therefore be interpreted as pure grav- 
itational waves. By separation of variables the equations governing these modes can 
be reduced to the 1-dimensional Schrodinger equation 

d 2 
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where the total potential Vi + V is dominated by Vi — 1(1 + l)r -2 , I = 2,3,... 
which is simply the effective potential for null geodesies with angular momentum 
L 2 = 1(1 + 1). This shows that the optical geometry is closely related to the trapped 
modes of gravitational radiation since the maxima and minima of Vi clearly correspond 
to the necks and bulges of the optical geometry respectively. We turn now to the 
analysis of the optical geometry of the double layer uniform density models. The 
equation of state is parametrized by the constant densities p+ and p_ of the interior 
and exterior zone respectively, as well as the transition pressure pt where the density 
is discontinous. Given these parameters, a stellar model is obtained by specifying 
the central pressure p c . Below we display the optical geometry and the axial mode 
potentials for three stars of this type. The parameter values are chosen to obtain three 
distinct examples of optical geometries with pronounced double necks, corresponding 
to two deep potential wells. As the wells are fairly rectangular-shaped with a certain 
width L and height Vq , one can make a crude estimate of the number of trapped modes 
by taking Z\ to be the wave function of an infinite square well, with the same width 
L, and count the number of states with energy less than Vq- The number of trapped 
modes is then approximately given by sfV^Ljix. 

The models considered here were chosen to have the same value of the tenuity 
a = R/M = 2.2749, and the overall scale was fixed by setting the total mass equal 
to unity (in units such that G a b — T a &). Below the optical geometries are displayed 
next to the potentials relevant for the axial modes with 1 = 2, the solid line indicating 
the total potential (V + Vi), the dashed line showing Vi and the dash-dotted line 
showing V. In the optical geometry the two shades of gray indicate the different 
density layers, and the white part corresponds to the exterior Schwarzschild geometry. 
The parameter values, given in units of the appropriate power of the total mass, are 
given in the respective captions. We may remark however that none of these models 
are likely to be stable, but that we have found models with less pronounced double 
necks that are probably stable. 




Fig. 1 The "peanut" geometry and the corresponding potentials. The parameter values are 
p+ = 2.512 x 10 3 , p- = 0.5024, p t = 35.89, and p c = 5.778 x 10 4 . 
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Fig. 2 The "bigup" geometry and the corresponding potential. The parameter values are 
p+ = 5.325 x 10 3 , p_ = 0.5047, p t = 29.58, and p c = 5.325 x 10 4 . 




Fig. 3 The "bigdown" geometry and the corresponding potential. The parameter values are 
p+ = 2.569 x 10 3 , p- = 0.5024, Pt = 36.69, and p c = 2.954 x 10 5 . 
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